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Abstract. The dynamics of randomly crosslinked liquids is addressed via 
a Rouse- and a Zimm-type model with crosslink statistics taken either from 
bond percolation or Erdos-Renyi random graphs. While the Rouse-type 
model isolates the effects of the random connectivity on the dynamics of 
molecular clusters, the Zimm-type model also accounts for hydrodynamic 
interactions on a preaveraged level. The incoherent intermediate scattering 
function is computed in thermal equiUbrium, its critical behaviour near the 
sol-gel transition is analysed and related to the scaling of cluster diffusion 
constants at the critical point. Second, non-equilibrium dynamics is studied 
by looking at stress relaxation in a simple shear flow. Anomalous stress 
relaxation and critical rheological properties are derived. Some of the 
results contradict long-standing scaling arguments, which are shown to be 
flawed by inconsistencies. 



1. Introduction 

Gelling liquids are part of everyday life. One encounters them, for example, when 
preparing a chocolate pudding or when sticking two materials together with the help of 
glue. From a microscopic point of view, gelling liquids consist of irregularly structured 
clusters of molecules or macromolecules. The formation of these clusters is either a result of 
intermolecular association, produced by e.g. van der Waals forces, electrostatic attractions or 
hydrogen bonding, or a result of chemical reactions such as polycondensation, polymerisation 
or vulcanisation induced by a chemical crosslinker [1,2]. Intermolecular association, also 
called physical gelation, leads to weakly bound clusters, which typically form and dissolve 
reversibly in the course of time during an experiment. On the other hand, chemical gelation 
leads to permanent clusters at temperatures of interest, and it is this situation that we will 
exclusively consider here. 

When increasing the concentration of crosslinks in a liquid (sol) one observes a more and 
more viscous behaviour under shear stresses, until a sudden transformation to an amorphous 
solid state takes place at a certain critical crosslink concentration. This point marks the 
gelation transition or sol-gel transition. The static shear viscosity diverges at the transition, 
and the onset of a static shear modulus is found. 
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Carothers [3] was the first to interpret the gelation transition as due to the formation of 
a macroscopic cluster of molecules in the system. His considerations were quantified and 
refined by Flory [4, 5] and Stockmayer [6, 7] to what is nowadays called "classical theory", a 
percolation model of tree-like structures, closely related to percolation on Bethe lattices [8]. 
So the classical theory arises [9] in the mean-field approximation of lattice-bond percolation 
[10]. Stauffer [11] and de Gennes [12] suggested the latter as a mathematical model for 
gelation, in particular, if caused by polycondensation. Lattice-bond-percolation clusters may 
also contain loops, and the spatial dimension becomes relevant, too. More importantly, upon 
identifying the gelation transition with the lattice-bond-percolation transition, it is revealed 
to be a continuous phase transition. Its driving parameter is crossUnk concentration, not 
temperature. Within this theoretical picture, the critical behaviour at the gelation transition 
is dictated by scaling and universahty [13, 10]. 

The resulting predictions for static properties of gelation clusters agree well with 
experiments in the vicinity of the sol-gel transition [14, 15] — a substantial improvement over 
the mean-field like classical theory. As far as dynamical phenomena are concerned, a variety 
of competing attempts have been made to seek an interpretation in terms of the percolation 
picture, see e.g. 1 16-18] for contradictory predictions concerning the shear viscosity. Yet, 
aU of these attempts rely on more or less ad hoc assumptions needed to compensate for the 
lack of thermal fluctuations or any sort of dynamics in a pure percolation model. Rather, 
the appropriate strategy should be to start from a (semi-) microscopic dynamical model 
for gelation clusters, from which the desired Unk to quantities in percolation theory can be 
deduced. This route will be followed here. Other analytical approaches to gelation from a 
microscopic model include [19-27]. Among others, they describe thermostatic fluctuations 
in the gel phase and calculate the static shear modulus. Computer simulations of microscopic 
models for gelation have been done by e.g. [28-34]. 

In this survey we wiU concentrate on the sol phase and report on results obtained in 
[35-42]. The dynamics of the sol phase is characterised by strong precursors of the gelation 
transition, even well below it. These include anomalous, stretched-exponential decays in time 
of both dynamical density correlations [43] and shear-stress relaxation [44]. Both decays are 
characterised by typical time scales which diverge when the critical crosslink concentration 
is approached. Our exact results on critical rheological properties contradict long-standing 
scaling arguments, which are shown to be flawed by inconsistencies. 

The paper is organised as follows. In Section 2 we briefly lay out a suitable generalisation 
of the usual Rouse and Zimm model for linear polymers to describe gelhng liquids. The model 
is then used to investigate time-dependent density fluctuations in Section 3. Section 4 deals 
with stress relaxation and critical rheological properties in a simple shear flow. Both Section 3 
and Section 4 are subdivided in a part pertaining to the Rouse model, a part pertaining to the 
Zinnm model and a part where the results are discussed and put in a wider perspective. FinaUy, 
Section 5 adds some closing remarks. 



2. Rouse and Zimm model for randomly crosslinked monomers 

In this section we give a brief description of a model which is to be considered a 
theoretical minimal model for the dynamics of gelling complex fluids. This model is a 
generalisation of one of the most fundamental models of polymer physics [45^8] to the 
case of randomly connected monomers. In this context, it has been discussed before by e.g. 
[49-56, 35, 37-39, 57-59, 40-42]. 
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2.1. Dynamical equation 

We consider N point-like monomers, which are characterised by their time-dependent 
position vectors Ri{t),i = 1, . . . , A'', in three-dimensional Euclidean space R^. Permanently 
formed, harmonic crossUnks connect M randomly chosen pairs of particles {ie,je)< where 
1 < ie je < N for all e = 1, . . . , M . The potential energy associated with these entropic 
Hookean springs takes the form 

e— 1 — 1 

where the length a > plays the role of an inverse crosshnk strength, and physical units 
have been chosen such that k^T = 1. It will be convenient to specify a given crosslink 
configuration Q := {{ie,je)}^i in terms of its N x TV-connectivity matrix F, which is 
defined by the right equality in (1). For part of what follows this setting could be generalised 
to the crosshnking of N identical molecular units which consist themselves of a given number 
of monomers that are connected in some fixed manner, such as N identical chains, rings or 
stars of monomers [37, 38]. For the ease of presentation, however, we will not consider such 
a generahsation here. 

We study the dynamics of these harmonically crosslinked monomers in the presence 
of an incompressible solvent fluid, which may induce hydrodynamic interactions between 
them. Hydrodynamic interactions will be incorporated on a preaveraged level in the spirit of 
Kirkwood and Riseman [60] and Zimm |46]. This is a traditionally accepted way of doing so 
albeit the limitations of this approach are still not sufficiently weU explored [47, 48]. We also 
aUow for the presence of an externally imposed, simple shear flow in a;-direction 

v(r,i) :=7(i)ye^ (2) 

with a time-dependent shear rate j{t). Here r = {x,y, z). A purely relaxational monomer 
dynamics is then described by [47, 48] 

- ^{R,(t), i) = - E H:,} + Ut) (3) 

for i = 1, . . . , n. This is the defining equation of the Zimm model for crosslinked monomers 
(in solution). The rest of this subsection is devoted to a brief explanation and discussion of 
(3), see [41,42] for more details. 

The jointly Gaussian thermal noises in (3) have zero mean and covariance 

£,i{t)i]{t') = 2H°'^ 5(t — t')l, as is required by the fluctuation-response theorem. As 
usual, the ^^ "thermalize" the system in the long-time limit. Here, the dagger denotes the 
transposition of a vector, 5 the Dirac-delta function and 1 the 3 x 3-unit matrix. 

Interactions between the monomers and the solvent fluid are subsumed in the spatially 
isotropic and homogeneous preaveraged mobiUty matrix 

^ Wj + {l-5ij)h{K^Tr/nij)]. (4) 



ueq 



It emerges [41,42] from taking Oseen's expression [61,60] for the mobility tensor and 
averaging it with respect to the suitably normalised Boltzmann weight ^ . However, 
when it is indispensable to have a positive definite mobility matrix in the sequel, we will 
replace the Oseen tensor with the Rotne-Prager-Yamakawa tensor [62, 63] in this procedure. 
Depending on which tensor is used, the function h in (4) is given by [64] 

h(x)-=l _ Oseen, 

\ erf (v^) — (1 — e ^)/0rx Rotne-Prager-Yamakawa. 
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The expression in the second line of (5) involves the error function erf and reduces to 

the expression of the Oseen case asymptotically as a; J, 0. The diagonal term in the 
preaveraged mobility matrix(4), which is proportional to the Kronecker symbol 5i^j, accounts 
for a frictional force with friction constant Q that acts when a monomer moves relative 
to the externally imposed flow field (2). The non-diagonal term reflects the solvent- 
mediated average influence of the motion of monomer j on monomer i. The parameter 
K := \/&JttC,/ (Q-K-qsa) involves the solvent viscosity rjs and serves as the coupUng constant 
of the hydrodynamic interaction. Formally setting k = in (4) yields H°'^ = C^'^Sij, and the 
Zimm model for crosshnked monomers reduces to the Rouse model for crosslinked monomers 
[35^0] 

I^Ri{t)-.{R,{t),t) = -\^^+m, (6) 

where i = 1, . . . , n and the jointly Gaussian thermal noises have zero mean and co variance 

^i{t)$}j{t') = {2/Q6{t — t')l. It is only for convenience that we introduced the Rouse 
model as the special case k = of the Zimm model here. Physically, it has its own standing 
as the minimal model for polymer melts under theta conditions, see e.g. [47, 48] for the case 
of linear polymer chains. In particular, all the approximations that entered the derivation of 
the (off-diagonal part of the) preaveraged mobility matrix H'^i do not affect the Rouse model, 
of course. 

It remains to explain the quantity TZij in (4), which is simply the mean squared 
displacement between monomers i and j in the thermal-equiUbrium state characterised by 
the suitably normalised Boltzmann weight ^ c^^ . In order to write down a formula for TZij, 
let us remark that, by construction, the connectivity matrix F = T{G) is block-diagonal with 
respect to the clusters of a given crossUnk configuration G (which are the maximal connected 
components of Q). Moreover, F(^) possesses as many zero eigenvalues as there are clusters in 
G. This is easily seen from the fact that the centre of mass of each cluster does not feel a force 
from the potential energy V. Hence, F cannot be inverted, but it possesses a Moore-Penrose 
pseudo-inverse Z [65], which is the inverse of F on the complement of its zero eigenspace and 
zero elsewhere. It can be represented as Z := (1 — Eo)/r, where Eq denotes the projector 
on the zero eigenspace of F in IR^ and 1 denotes the N x A/^-unit matrix. The mean-squared 
displacement TZij is then given in terms of Z according to 



^ f Zi^i + Zjj — 2Zij if i and j belong to the same cluster, 

' \ -|-oo otherwise. 

There is also another interpretation for TZij , which we will use below: Viewing each monomer 
as an electric contact and each crosslink as a unit Ohmian resistor connecting two contacts, 
TZi,j is the effective electric resistance between the contacts i and j of this corresponding 
electrical resistor network [66]. This exact correspondence between Hookean bead-spring 
clusters and Ohmian electrical resistor networks relies on the Unearity of Hooke's and Ohm's 
law. 

Since both the connectivity matrix F and the preaveraged mobility matrix H^'^ are block- 
diagonal, it follows that clusters move independently of each other in this model. The saUent 
feature of the Zimm and Rouse equations (3) and (6) is that they are linear in the monomers' 
positions. Hence, they admit an exphcitly known solution. The results we present in this 
paper rely heavily on this solution. 
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2.2. Average over crosslink ensemble 

So far, everything in this section was meant for an arbitrary but fixed realisation Q of 
M crosslinks among N monomers. For practical reasons, Q can never be determined 
experimentally in macroscopicaUy large gelUng fluids. Neither should physicaUy meaningful 
observables depend on specific microscopic details of G, but only on some macroscopic 
characteristics of it. Therefore, we follow the general philosophy of the theory of disordered 
systems and take ^ as an element of a statistical ensemble of crosslink configurations, within 
which it occurs with probability Pn{Q)- The just made statement on physically meaningful 
observables A(Q) now translates into a self-averaging property: the two quantities A(Q) and 
its ensemble average Pn{Q')-A{Q') coincide (with probabihty one) in the macroscopic 
Umit. Therefore we will compute the macroscopic Umit 

{A) := lim y^PN{g)A{g) (8) 

AT— >oo ^ — ' 
Q 

of such averages with a fixed crosslink concentration c := limjv->oo M/N . This will be done 
for two different crosshnk ensembles. 

(i) Clusters are generated according to three-dimensional continuum percolation, which 
is closely related to the intuitive picture of gelation, where monomers are more likely to 
be crosslinked when they are close to each other. Since continuum percolation and lattice 
percolation are beUeved to be in the same universality class [10], we employ the scaUng 
description of the latter. It predicts [10] a cluster-size distribution of the form 

T„ ^ exp{— n/n*} (9) 

for £ := (Ccrit — c) <C 1 and n ^ oo with a typical cluster size n*{e) ^ e"^/'^ that diverges 
as e ^ 0. Here, a and r are (static) critical exponents, see Table 1 below for their numerical 
values. 

(ii) Each pair of monomers is chosen independently with equal probability c/N, 
corresponding to Erdos-Renyi random graphs, which are known to resemble the critical 
properties of mean-field percolation [9]. After performing the macroscopic Umit, there is no 
macroscopic cluster for c < Ccrit = 1/2 and almost all clusters are trees [67]. Furthermore, 
all n"~^ trees of a given "size" n, that is, with n monomers, are equally likely. The cluster- 
size distribution can also be cast into the scaling form (9) with the exactly known critical 
exponents r and a listed below in Table 1. 

3. Time-dependent density fluctuations 

In this section we address dynamical properties of gelhng liquids in thermal equihbrium. 
Therefore we will assume throughout this section that there is no externally imposed shear 
flow, i.e. 7 = 0. 

Experiments [43, 68] on quasi-elastic light scattering in gelling liquids allow to measure 
how spatial density fluctuations of a given wave vector q are correlated to each other at 
different times t. This information is encoded in the incoherent intermediate scattering 
function 



S(q,t):= lim — e'9-[Ri(«+to)-Hi(to)] . (lO) 

2—1 

The right-hand side of (10) is determined by the solution Ri^t) of the dynamical equation (3) 
for a given crosslink reaUsation Q and with initial conditions being imposed at time Iq. The 
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average over the thermal noise and the subsequent limit to — > — c» in (10) ensure that the 
system reaches its thermal-equiHbrium state. Then, for large retardation times t, one expects 
[69, 43] that this correlation is determined by the slowest relaxation processes in the system. 
Due to the independent motion of different clusters in the model under consideration, the 
slowest relaxation processes correspond to the centre-of-mass diffusion of whole clusters of 
monomers. This argument can be quantified — see e.g. [35], [41] or Eq. (4.12) in [37] — and 
yields 

S{q,t) e^^{-qHD{Nu)} . (11) 

fc=i 

Here we have set (j : = | q | and introduced the clusters Mk ,k = 1 , . . . , if , of the given crossUnk 
configuration Q. The number of monomers in the cluster Mk is denoted by Nk and 



(12) 



defines its diffusion constant in terms of the mean-square displacement of its centre of mass 
-RcM fc {t) := J2ieXk '^Sht equality in (12) follows from a short calculation 

with the exact solution of the dynamical equation (3). It was previously established in [70]. 
Another diffusion constant has been introduced by Kirkwood [48, 47] 

It provides an upper bound to the former, 

D{Uk) < D{Nk) , (14) 

as can be shown by applying the Jensen-Peierls inequality, see e.g. Sect. 8c in [71], to (12). 
Customarily, one also defines an effective diffusion constant i'eff for the whole gelhng liquid 
by 

^i'-iL"i»7„ <"M.,.)-gj^^. (15) 

Since S{q, t) is expected to develop a time-persistent part in the gel phase, Des is expected 
to vanish when approaching the gelation transition from the sol side. 

3.1. Rouse dynamics 

We recall from Sect. 2.1 that in the absence of hydrodynamic interactions, k = 0, we have 



constant (13) are equal 



^ ^ ^hj- Hence, the cluster-diffusion constant (12) and the Kirkwood diffusion 



D(A4) = 5(A4) = 7^ , (16) 

and inversely proportional to the number of monomers in the cluster [35]. In other words, 
cluster topology does not influence diffusion within Rouse dynamics. 

Next, we discuss the long-time behaviour of the incoherent intermediate scattering 
function in the macroscopic hmit. According to Sect. 2.2, this amounts to calculating the 
average of (11) 

{S{q, t)) (E ^ eM-Q'tD{Afk)}) . (17) 
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Thanks to (16) this average is easily performed by reordering the clusters according to their 
size 



(5(q,i))*~~^nr„e-«^*/(«'^), (18) 

71=1 



where 



fe=i 

is the cluster-size distribution and (18) holds in the absence of an infinite cluster. Using the 
scaling form (9) of r„, we find [35, 37J 

(^-|^^ { [Vi*(£)](^-i/2)/2exp{-const. e>Q. ^^^^ 

At the critical point, the long-time decay is algebraic with a critical exponent y = r — 2. In 
the sol phase one has a Kohlrausch or stretched-exponential behaviour with a time scale that 
diverges as tq{e) ^ (C/?^)^ '' with a critical exponent = 1/cr, when the critical point is 
approached. 

For the effective diffusion constant (15) we conclude from (20) that it vanishes like 

(£>eff) ^ \ira\q\{e)\ "^^"^^ ~ e" with a = (3 - t)/g (21) 

as £ i 0. 

The exponent a could have also been deduced directly from the right expression in (15). 
Indeed, given any cluster-additive observable, a reordering of the clusters according to their 
size yields 

= (E ^ ^^^^)) = E "^"(^)" ' (22) 



where 



1/^1 \ 
(^)„:=-(5]-5;v„n^(^fe)) (23) 



N 
fe=i 



is the partial average of A over all clusters of a given size n. Now, if the partial averages 
exhibit the critical divergence 



then 



{A) ~ as e i with u = {2-t-\- b)/a, 
provided that u> 0. 



(25) 
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Table 1. Numerical values for the critical exponents of the cluster-size distribution (9) and the 
two fractal dimensions of Gaussian phantom clusters in (27). The values are listed for cluster 
statistics according to three-dimensional bond percolation (3D) and Erdos-Renyi random 
graphs (ER). 





T 


a 


ds 




3D 
ER 


2.18 
5/2 


0.45 
1/2 


1.33 
4/3 


3.97 
4 



3.2. Zimm dynamics 

In contrast to the free-draining limit described by Rouse dynamics in the last subsection, 
one expects that with hydrodynamic interactions being present, cluster topology will have an 
influence on the diffusion constants. 

For simplicity, let us start with the Kirkwood diffusion constant. In order to extract a 
size dependence out of D, we look at the average over all clusters of a given size n 
and study its behaviour as a function of n. More specifically, we will perform this average 
precisely at the critical concentration Ccrit, where we expect an algebraic decrease as n oo 
due to the absence of any other length scale at criticaUty. Indeed, using the Oseen tensor for 
the hydrodynamic interactions we deduce from (13), (4) and (5) that 

where A is some dimensionless proportionality constant. The asymptotic behaviour of the 
average over the resistances in (26) is derived in [41]. The derivation has to distinguish 
between the two different cases for the crosslink ensemble. For Erdos-Renyi random graphs 
the asymptotics can be deduced from the exact probabiUty distribution of TZi^j in [72]. For 
three-dimensional bond percolation we use the scaling form of the probability distribution, 
which was established within two-loop order of a renormalisation-group treatment of an 
associated field theory [73, 74]. Equation (26) involves the fractal Hausdorff dimension 

df^ := 24/(2 - ds) (27) 

of Gaussian phantom clusters, which also determines the scaUng of their radius of gyration 
according to [50, 52, 55] 

^ y URi-Rjf) I 

2j^2 / ^ \ ^ ' J ' I n\ c=Cci 



RgyT,n 



f . (28) 



The other fractal dimension in (27) is the spectral dimension ds of the incipient percolating 
cluster [75,76]. Their numerical values are listed in Table 1. We conclude from (26) that 

Dn shows a crossover from Rouse behaviour £)„ for n < n(/t) K~^^^^~^^'^f ^ to 

Zimm behaviour 



Dn^n-'/^'f' -l/i?gyr,„ (29) 



for asymptotically large n > n{K). 

Now we turn to the averaged diffusion constant 

£>„:=(D)„|^_^ "~°°n-^° (30) 

' C — Cent 
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of clusters of size n at the gel point, which is also expected to obey a critical scaling for large 
cluster sizes n. From the Jensen-Peierls inequality Z)„ < Z)„, see (14), we then infer the 
inequality 

bo > (31) 

for the critical exponents. Figure 1 shows numerical data for the cluster diffusion constant 
Dn, plotted against n, for different values of the hydrodynamic interaction strength k. The 
crosslink ensemble in Fig. 1(a) corresponds to Erdos-Renyi random graphs. In Fig. 1(b) 
crosslinks were chosen according to three-dimensional bond-percolation. In the numerical 
computations we have used H^^ corresponding to the Rotne-Prager-Yamakawa tensor so 
that a positive definite mobility matrix is always guaranteed. Like the Kirkwood diffusion 
constant, £)„ also exhibits a crossover from Rouse to Zimm behaviour at a cluster size 
comparable to n(K). Figure 2 shows the exponent bo of the power-law fit (30) to the data 
of Fig. 1 in the large n-regime for the different values of k. The horizontal dashed lines 
in Figs. 2(a) and (b) correspond to the exponent value of the respective Kirkwood 

diffusion constant. The bigger exponent values that occur for small values of k still show 
residual Rouse behaviour for the largest system sizes we treated. For bigger values of k the 
crossover can hardly be felt any more in the largest systems, and the extracted exponent value 
bo corresponds to Zimm dynamics. This value is very close to the scaUng exponent in (29) 
for the Kirkwood diffusion constant, and, in fact, we conjecture that 

bo = l/df\ (32) 

We now turn to the long-time behaviour of the incoherent intermediate scattering 
function (10). The asymptotics (11), (22) and Jensen's inequality yield the lower bound [41] 

oo 

>5^nr„e-«'*^''. (33) 

n=l 

In fact, there is numerical evidence that this inequality actually captures the correct long-time 
asymptotics of {S{q,t)). Evaluating the right-hand side of (33) for large times t, this then 
leads to the scaling form [41]| 

{S{q,t)) ~ (^-^^ I [i/i*(£)].te-i/2)g^p|_^oj^gt_[^/^*(g)j.| (34) 

with the time scale t*{£) g~^e~^. The exponents are given by 

x = {l + bD)-\ y={T-2)/bD, z = bD/a (35) 

and are expressed in terms of bo ^ 0.25, see (32) and Table 1. The Rouse Umit (20) of (34) 
corresponds to setting bo = 1 in the above expressions. 
The critical vanishing 

(£>eff ) ~ with a= {2-T + bD)/o- (36) 

of the effective diffusion constant follows from directly from (22) - (25) provided that a > 0. 
This condition is fulfilledfor three-dimensional bond percolation where a w 0.16, but violated 
for Erdos-Renyi random graphs. Finally, we like to point out that, regardless of the cluster 
statistics, the ensemble averaged diffusion constant {D) never vanishes at the critical point. 
This is simply because it has non-vanishing contributions from all clusters, which add up. 

i Note that there is a misprint in the second line after Eq. (30) in [41]. The algebraic prefactor in the scaling form of 
the function s(A) should read \^(y-'^/'^'> instead of 
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Figure 1. (a) D„ at the gel point for mean field percolation and different hydrodynamic 
interaction strengths, (b) Same for three-dimensional bond percolation. 
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Figure 2. (a) Critical exponents ho, corresponding to a power-law fit Dn ~ n"**" in 
Fig. 1(a). (b) Same for Fig. 1(b). 



3.3. Discussion 

We have studied the critical scaHng _D„ ^ n^''° of the averaged cluster diffusion constants 
over clusters of size n and used it to obtain the scaling behaviour of the intermediate 
incoherent scattering function {S{q, t)) near criticaUty. The associated critical exponents 
are summarised in Table 2. Within Rouse dynamics cluster diffusion constants are inversely 
proportional to the cluster size n, irrespective of the cluster topology, that is, 6d = 1. Zimm 
dynamics leads to bo = see (32), and topology does play a role: Indeed, it is 

well known [48] that within Zimm dynamics the diffusion constant of a linear chain of n 
monomers decreases as n~^/^. Since bn ~ 0.25 < 1/2, this means that, on average, a 
monomer in a branched cluster feels less friction — ^which is intuitively appealing, because 
monomers in the interior of a cluster should be dragged along. Second, (28), (30) and (32) 
imply for Zimm dynamics that Dn ~ l/i?gyr,n. Hence, this relation does not only hold for 
linear chains, for which it has been well known [48], but in an average sense for all percolation 
clusters. 

Concerning the scaling exponents of the incoherent intermediate scattering function. 
Table 2 shows that neither Rouse nor Zimm dynamics provides even a reasonably good 
description of the experimental findings, despite their strong scatter There are several reasons 
for the discrepancies between the model predictions and experiments, (i) Our results 
pertain to ^-conditions, in so far as excluded- volume interactions have been neglected in the 
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Table 2. Summary of critical exponents for cluster diffusion constants and the incoherent 
intermediate scattering function (see Eqs. (30), (20) and (34) for their definitions). The 
numerical values for Rouse and Zinun dynamics — Usted for cluster statistics from three- 
dimensional bond percolation (3D) and Erd6s-R6nyi random graphs (ER) — are compared to 
experimental findings. 





Zimm 


Rouse 










Exponent 


3D 


ER 


3D 


ER 


[68] 


[77] 




[78] 


ho 


0.25 


1/4 


1 


1 










X 


0.80 


4/5 


1/2 


1/2 


0.66 


0.3- 


0.8 


0.64 


y 


0.71 


2 


0.18 


1/2 


0.27 


0.2- 


0.3 


0.34 


z 


0.56 


1/2 


2.22 


2 


2.5 








a 


0.16 


(*) 


1.82 


1 


1.9 


0.5- 


1 


1.9 



(*) no divergence 



models. Excluded-volume interactions could cause a swelling of the clusters, which results 
in a different fractal Hausdorff dimension, (ii) We chose cluster statistics according to 
three-dimensional bond percolation. This accounts well for crosslinldng in a dense melt, say, 
but not in dilute solutions, (iii) It has been suggested [68] that hydrodynamic interactions 
between monomers in a cluster are screened by smaller clusters in the reaction bath so that 
the Rouse rather than the Zimm model should apply. Our analysis supports this conclusion 
in so far as the exponents of the Rouse model are closer to the experimental values. So the 
more striking failure of the Zimm model can be traced back to a too slow decay of D„ with n. 
(iv) Preaveraging of the hydrodynamic interactions is an uncontrolled approximation, and it 
remains to be seen what a full treatment of hydrodynamic interactions predicts for the critical 
dynamics of gelling solutions. 

4. Stress relaxation 

Gelling liquids exhibit striking rheological properties which have been continuously studied 
over the years by experiments [79-86], theories [11,16, 50, 87, 88, 35, 89, 36, 37, 90, 38, 39] 
and simulations [28,29, 31,91, 32-34]. For example, when subjected to the homogeneous 
shear flow (2), distinct relaxation patterns are observed, which are due to the participation of 
many different excitation modes of all sorts of clusters. IVIore precisely, experiments suggest 
the scaling form [92, 44, 80-82, 86, 87] 

{G{t)) - fr^g{t/t) with t{e) ~ (37) 

for the macroscopic (shear-) sttess-relaxation function in the sol phase for asymptotically 
long times t and crosslink concentrations close to the critical point, i.e. for e <C 1. The typical 
relaxation time t diverges with a critical exponent z > for e J, 0. The scaling function g is 
of order unity for small arguments so that one finds the algebraic decay {G{t)) ~ t~'^ with a 
critical exponent < A < 1 for f ^ oo at the critical point. For large arguments, g decreases 
faster than any inverse power. Sometimes a stretched exponential has been proposed for g in 
this asymptotic regime [82, 87]. 

In this section we will investigate to what extent such critical properties can be predicted 
by the Rouse and the Zimm model. Thus we will explore the consequences of the dynamics 
(3), resp. (6), in the presence of the externally appUed simple shear flow (2). In reaction to 
the flow, the system of crosslinked monomers builds up an intrinsic shear stress. FoUowing 
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Kirkwood, see e.g. Chap. 3 in [48] or Chap. 16.3 in [47], this shear stress is given by the force 
per unit area exerted by the monomers 



N 



N 



(38) 



i=l 



Here, Hi (t) is the solution of the equation of motion (3) with some initial condition at time 
^0 in the distant past (so that the noise average yields a thermahzed state in which all transient 
effects stemming from the initial condition have died out). Moreover, po stands for the 
monomer concentration and Fi{t) := — dV/ dRi (t) is the net spring force acting on monomer 
i at time t. The explicit computation [37, 38] of the right-hand side of (38) yields 



(T{t) = G(0) 1 + J dt' G(t - t') 7(t') ^ 



2 //,ds7(s) 1 








(39) 



for arbitrary strengths of the shear rate 7(f). Here, we have defined the stress-relaxation 
function 



G{t) := ^ Tr 



(1 - Eo)exp 



(40) 



as a trace over the matrix exponential of f := (H®i)^/^r (H'^i)^/^. Due to the occurrence of 
the spectral projector Eq on the kernel of F, this trace is effectively restricted to the subspace 
of non-zero eigenvalues. 

For a time-independent shear rate 7, the shear stress (39) is also independent of time. 
The viscosity r} is then related to shear stress via 



7P0 Po io 3 2iV 



1-En 



(41) 



Apparentiy, the viscosity is determined by the trace of the Moore-Penrose inverse of F. The 
normal stress coefficients are given by 



^(1) ,-_ 



7^Po 



Y2 



and 



^(2) ._ ^VjV. 



= 0. 



(42) 



(43) 



The vanishing of \E''^' is typical for Rouse/Zimm-type models and has been well known for 
the case of Unear polymers [48]. Since F is block-diagonal with respect to the clusters, the 
observables G{t), rj and ^'^^^ are all cluster-additive in the sense of (22). 

The scaling form (37) of the macroscopic stress-relaxation function {G{t)) implies that 
the macroscopic viscosity and first normal stress coefficient exhibit a critical divergence 

(??) ~ and (*(^)) ~ (44) 

at the sol-gel transition as £ J, with critical exponents given by the scaling relations [44, 39] 

= 2(1 - A) and t = z{2 - /^) = k + z . (45) 

Thus, it suffices to know any two of the four critical exponents ^,z,k and L 
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4.1. Rouse dynamics 

For Rouse dynamics we have V = T/( so that the computation of the stress-relaxation 
function, the viscosity or the first normal stress coefficient requires the knowledge of spectral 
properties of the cormectivity matrix F. 

Concerning the macroscopic viscosity (77), there are several ways of calculating the 
critical exponent k in (44). The different ways explore connections to problems in different 
branches of research. Given a cluster A4, the trace of the Moore-Penrose inverse of r(A4) 
can be expressed in terms of the resistances (7) according to [36, 37] 

-Eo(7^fc)" 



'^i'j- (46) 



We stress that this is an exact relation [66]. It has nothing to do with electrical analogues 
put forward in scaling arguments [69]. For the case of Erdos-Renyi random graphs there are 
only tree clusters for c < Ccrit = 1/2. In this special case the resistance TZij reduces to the 
graph distance of i and j in Afk, and the right-hand side of (46) is known as the Wiener index 
WiMk) in graph theory. From a graph-theoretical point of view, the right equality in (46) 
follows also as an application of the matrix-tree theorem, see e.g. [93], Thm. 5.5. Moreover, 
the partial averages {W) are exactly known [72], and, using (22), one finds the exact result 
[36,37] 



24c 



(47) 



It can be interpreted as a critical divergence with exponent k = 0. Alternatively, (47) can also 
be obtained from a replica approach [37] instead of using graph theory. The replica approach 
is also capable of providing us with higher inverse moments {N~^ Tr [(1 — Eq) /T"]) for not 
too large positive integers v [90]. Using these results for i/ = 2, a (somewhat lengthy) exact 
expression for ^P^^^ was derived in [39] for crosslink statistics from Erdos-Renyi random 
graphs. It exhibits the critical behaviour 

(4f(i))^e-^ with i = 3. (48) 

Now we turn to the crosslink ensemble of three-dimensional bond percolation. In order 
to proceed from (46) in this case, one needs to know the average resistance {'R-ij)n between 
two nodes in bond-percolation clusters of size n. Luckily, random electric resistance networks 
have been studied extensively, and the asymptotic behaviour 

(7^^J)„ ~ n^" with := (2/4) - 1 (49) 

can be extracted [36,37] from highly developed renormalisation-group treatments of an 
associated field theory [73,74]. Thus, (46), (22) and (25) lead to the critical behaviour 
(77) ~ with 

k = {l-T + 2/d,)/a (50) 

as £ J, 0. Of course, this exact scaling behaviour reduces to the Erdos-Renyi result A; = 
from (47), when inserting the appropriate mean-field values for the exponents. 

None of the above approaches is able to yield any of the other critical exponents A and 
z — or also £ in the case of three-dimensional percolation statistics. Here, a connection to 
random walks in random environments is helpful. For the time being, let us concentrate on 
the case of three-dimensional percolation statistics, where the maximum number of bonds 
emanating from any vertex is hmited to m = 6 on the simple cubic lattice. Now, consider a 
random walker — coined "blind ant" by de Geimes [12] — that moves along a bond from one 
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site to another in the same cluster at discrete time steps [10, 76, 94,95]. If the ant happens to 

visit site i at time s, which is connected with m,; < m bonds to other sites, then it will move 
with equal probabihty 1 /m along any one of the mj bonds within the next time step and stay 
at site i with probabihty 1 — mj /m. By definition of the connectivity matrix T of the cluster, 
one has Ta = m.i for its diagonal matrix elements, = — 1 if two different sites i ^ j are 
connected by a bond and zero otherwise. Hence, the associated master equation for the ant's 
sojourn probabihty Pi{s) for site i at time s reads 

p,{s + 1) = (1 - ru/m)p,{s) + ^(-ri,/m)p,(s) , (51) 

which is equivalent to 

Pi{s + 1) - pi{s) = ^ Tijpjis) . (52) 

j 

Here the summation extends over all sites in the cluster. For long times s 3> 1, it is legitimate 
to replace the difference (quotient) on the left-hand side of (52) by a derivative. This yields 
the solution p,;(.s) = [c^"'"/™] which corresponds to the initial condition pi(0) = Si^i^. 

Next we consider (s) := {pi^ (s))„ |^^q, the mean return probabihty to the starting point 
after time s, where the average is taken over all critical percolation clusters with n sites. 
Clearly, these definitions are independent of the starting point io, because on average there is 
no distinguished site by assumption. Thus we can also write 



p(")(s) = (1 Tre-**^/" 



(53) 

=0 

for finite n. The return probabihty behaves as [76, 95, 94] 

(s) ~ s-'^'/^J^{s/s„) + 1/n , (54) 

where s„ ~ n'^^'^' and the cut-off function J^{x) is of order one for x < 1 and decreases 
rapidly to zero for x ^ oo. Basically, (54) says that for times s ^ s„ the walker has no 
memory of where he had started from. For times s < Sn the fractal-like nature of a cluster at 
c = Ccrit leads to an algebraic decrease of the return probability, which involves the spectral 
dimension ds. Now, assuming that {G{t)) obeys the scaling form (37), the information 
provided by (54) for c = Ccrit is sufficient to conclude [40] the exponent relations 

A=^(r-1) and z = . (55) 

When plugging (55) into (45), we recover (50) and get the new scahng relation 

£={l-T + 4:/ds)/(7. (56) 

Since, the critical behaviour of Erdos-Renyi random graphs coincides with that of mean-field 
percolation, we get the missing exponents A and z for that case by inserting the mean-field 
values into (55). 



4.2. Zimm dynamics 

The matrix F, which determines stress relaxation, is by far more complicated than F in the 
presence of hydrodynamic interactions. In particular, it reflects cluster topology only in a 
much more subtle way than F. In fact, it was that apparent encoding of topology in F that 
made the analytical methods of the last subsection work. In the absence of suitable analytical 
tools, numerical methods remain to investigate stress relaxation in the Zimm model. 
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Figure 3. Numerical data to determine the scaling (57) for random clusters in the case 
of Erdos-Renyi random graphs (left column) and three-dimensional bond percolation (right 
column). In each case the averaged viscosity r/n (top) and normal stress coefficient "jji^^ 
(middle) are plotted for different strengths of the hydrodynamic interaction pai'ameter k as a 
function of the cluster size n on a double logarithmic scale. Power-lavt' fits to the data yield 
the exponents and 6^ as a function of k (bottom). 



We determined the scaling as n 



\e=0 



oo of the partial averages 

and 'fiP 



le=0 



(57) 



at criticality by numerically diagonalising F and performing the disorder average over the 
crosslink ensemble [42]. The critical exponents k and £ then follow from (25). All numerical 
computations were done with the Rotne-Prager-Yamakawa tensor for the hydrodynamic 
interactions. The reader who is interested in more details of the numerical computations 
is referred to [42]. 

In Figs. 3(a) and (b) we plot rj^ and '^n^ as a function of n on a double-logarithmic 
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scale for different values of the hydrodynamic interaction parameter k. Crosslink statistics 
were chosen according to Erdos-Renyi random graphs. The exponents 6^ and 6^ are obtained 
from power-law fits in the large n-range and are displayed in Fig. 3(c). The viscosity exponent 
decreases from 6^ = 0.28 for k = 0.05 to 6,, = 0.11 for k = 0.3. We recall from (46) and 
(49) that the Rouse exponent for k is exactly given by 6,, — 1/2. The exponent 65, of the 
normal stress coefficient ranges from 6* = 1.2 for k = 0.05 to 6* = 0.73 for k = 0.25. The 
exact Rouse value 6* = (4/ds) — 1 = 2 for k = follows from (48) and (25). 

The same is done for three-dimensional bond percolation in the right column of Fig. 3. 
Figures 3(d) and (e) contain rjn and \ respectively, as a function of n on a double- 
logarithmic scale for different values of k. The exponents br, and 6*, extracted by fitting 
the curves in Figs. 3(d) and (e) to a power law for large n, are shown in Fig. 3(f). The 
numerical values for 6^ are nearly identical to those obtained for Erdos-Renyi random graphs. 
Again, one observes a decrease from 6^ = 0.21 for k = 0.05 to 6^ = 0.11 for k = 0.3. 
The exponent 65, of the normal stress coefficient ranges from &^ = 1.1 for k = 0.05 to 
bq, = 0.78 for K = 0.25. The corresponding exact Rouse values 6^ = {2/ds) — 1 « 1/2 and 
bsij — (4/(is) — 1 ~ 2 for K = follow from (50) and (56) in the last subsection together with 
(25). 

A careful analysis of the data in [42] reveals that the true Zimm exponents 6^ and b^ 
are universal in k, and that their seeming dependence on k in Figs. 3(c) and (f) is most likely 
due to finite-size effects. More precisely, for small k the data suffer from a crossover to their 
respective Rouse values so that they come out too large. For large k, on the other hand, the 
asymptotics h{x) ~ 1 — (ttx)"^/^ as a; —> 00 of the lower line in (5) leads to a slower growth 
of brj and at intermediate n. Hence, the exponents come out too small for larger k. The 
most reliable values for the universal Zimm exponents br, and 6* should be obtained from 
around k w 0.3. It it these values which are listed in Table 3 below. The critical behaviour 
of the averaged viscosity {1]) ~ e~'' and of the averaged first normal stress coefficient 
(vl/(i)^ ^ for a polydisperse gelling solution of crosslinked monomers then follows from 
(25). For the viscosity this implies a finite value at the gel point for both, Erdos-Renyi random 
graphs and three-dimensional bond percolation. In contrast, the first normal stress coefficient 
is found to diverge with an exponent that depends on the cluster statistics. Choosing the 
cluster statistics according to Erdos-Renyi random graphs, we find £ « 0.54. The case of 
three-dimensional bond percolation leads to the higher value £ « 1.3. These exponent values 
are less than a third in magnitude than the corresponding exact analytical predictions of the 
Rouse model from (56) with the corresponding cluster statistics. AH exponent values are 
summarised in Table 3. 

4.3. Discussion 

A fairly complete scaling picture of the gelation transition has been obtained within Rouse 
dynamics. All critical exponents k, I, A and z of the stress-relaxation function in the sol phase 
and at criticality could be expressed in terms of two independent static percolation exponents 
a and r plus the spectral dimension of the incipient percolating cluster, see the scaling 
relations (50), (55) and (56). These scaling relations and the resulting numerical exponent 
values listed in Table 3 contradict the predictions k = 2^ — /3 and A = dv/ldv + k) of earlier 
scaUng arguments [80,81, 84, 16, 87, 104]. What is the reason for this discrepancy? The 
scaling arguments involve the fractal Hausdorff dimension df := d — P/v of rigid percolation 
clusters at Ccrit- Rouse clusters, however, are thermally stabilised, Gaussian phantom clusters 
with the fractal Hausdorff dimension d^^\ see (27) [50, 52, 55]. The latter is different from 
rf/ in space dimensions below the upper critical dimension rf„ = 6. Indeed, if one replaces df 
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Table 3. Summary of critical exponents for stress relaxation (see Eqs. (37), (44) and (57) for 
their definitions). The numerical values for Rouse dynamics are based on the scaling relations 
(49), (50), (55) and (56). Those for Zimm dynamics are based on the data analysis of Fig. 3. 
The values are hsted for cluster statistics according to three-dimensional bond percolation (3D) 
and Erd6s-R6nyi random graphs (ER), and are compared to some experimental findings. 

Zimm Rouse 

Exponent 3D ER 3D ER [96] [83] [97] [98] [99] [100] [101] [81] [84] [102] [103] 

k (*) (*) 0.71 0(#) 0.2 0.7 0.82 1.1 1.27 1.3 1.36 1.4 >1.4 6.1 

I 1.3 0.54 4.1 3 

A 0.79 1 0.72 0.71 0.69 0.67-0.68 0.66 0.70 0.33 0.69-0.77 

z 3.3 3 2.9 2.67 

br, 0.11 0.11 0.50 1/2 

6* 0.77 0.77 2.0 2 

no divergence logarithmic divergence 

(G) 

by dj in these scaling arguments, as one should consistently do within a Rouse description, 
the results will coincide with the ones obtained here. 

Since the long-standing scaUng relations k = 2i/ — and A = dv/ (du + k) involve the 
Hausdorff fractal dimension dj of rigid percolation clusters, it is sometimes argued that they 
describe the behaviour of a more realistic model, which, in addition to the interactions of the 
Rouse model, accounts for excluded-volume effects, too, see e.g. [104]. As far as we know, 
this claim has not been verified by analytical arguments within a microscopic model. One 
may even have doubts whether this claim is generally true: Extensive molecular-dynamics 
simulations [29] of a system of crossUnked soft spheres in three dimensions, with cluster 
statistics from percolation and an additional strongly repulsive interaction at short distances, 
yield the values k ~ 0.7 and A « 0.75, which are remarkably close to the predictions 
of the Rouse model for randomly crossUnked monomers, see Table 3. On the other hand, 
simulations of the bond-fluctuation model in [28] imply k Ri 1.3 and are thus in favour of the 
claim. However, the viscosity is not measured directly in these latter simulations. Rather it is 
derived from the scaling of diffusion constants and an additional scaling assumption that may 
be questioned [29]. Hence, it is an open problem to what extent the critical Rouse exponents 
of Table 3 are modified by excluded-volume interactions. 

In the context of dynamical critical phenomena, one usually expects dynamical scaUng 
to hold. Thereby one can infer critical properties of the gel phase from those of the sol 
phase. In particular, the critical behaviour of the shear modulus Go ~ |£|'^ follows from the 
scaling form (37) of the stress-relaxation function. The result /x = Az = (r — l)/fT involves 
only the two exponents cr and r of the cluster-size distribution. Using well-known scaling 
relations of percolation theory, this can be rewritten as ii = dv in terms of the correlation- 
length exponent v and the spatial dimension d. It is in agreement with the simple scaUng 
argument based on dimensional analysis of the free-energy density. In a recent letter [105], 
the scaling of entropic shear rigidity was analysed for both phantom chains and those with 
excluded-volume interactions. In both cases the gel was prepared by crosslinking a melt 
of chains with excluded-volume interactions. Our choice of percolation statistics combined 
with Rouse dynamics should be comparable to phantom chains prepared in an ensemble with 
excluded-volume interactions. However, the results of [105] for /it disagree with the above 
dynamic-scaling argument. The reasons for the discrepancy are not understood. 

Let us return to the sol phase and discuss the Zimm results, which are based on the 
exact numerical determination of the scaling exponents &^ and b^^, for the fixed-size averages 
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(57) of the viscosity and of the first normal stress coefficient. The resulting finiteness of the 
macroscopic viscosity (77) at the transition is clearly the most serious drawback of the Zimm 
model for randomly crosslinked monomers. Also, this failure comes unexpected, because a 
well-known scaling argument [14, 106, 50] predicts a logarithmic divergence. This scaling 
argument uses the (correct) scaling Z?„ ^1/ Rgyr.n of diffusion constants together with the 
Stokes-Einstein relation and yields 6^ = d/df — 1. Consequently, one gets from (25) the 
scaUng relation k = (1 — r + d/df)/a. Inserting hyperscaling and the fractal dimension 
of rigid percolation clusters, one would get fc = from that, which was interpreted as a 
logarithmic divergence. But as we remarked already earlier on in this subsection, the correct 
fractal dimension df for Gaussian phantom clusters is (ZyP\ For both cluster statistics this 
would give an unphysical negative value around — 1/4 for which can be definitely ruled 
out by our data.§ Thus, we conclude that the scaUng approach of [14, 106,50] does not 
apply to the Zinrni model for randomly crosslinked monomers. Another scaling approach 
to this model by [17] is also falsified by our data. On the other hand, Brownian-dynamics 
simulations of hyperbranched polymers were performed in [91]. They also account for 
fluctuating hydrodynamic interactions corresponding to k = 0.35, as well as for excluded- 
volume interactions and lead to hr^ = 0.13. This result is remarkably close to our finding 
hr, « 0.11 for the highest coupling strength k = 0.3 that we have considered, whereas 
experimental findings (see below) are consistently above our value. 

Next, we comment on how the Rouse and Zimm predictions for stress relaxation compare 
to experimental reality. Table 3 shows an enormous scatter of the experimental data. Thus, a 
serious check of theoretical predictions is currently severely hampered. The origin of this wide 
spread of the data is unclear so that even the question arose, whether the dynamical critical 
behaviour at the gelation transition was indeed universal [107]. Possible explanations for non- 
universal behaviour include the spUtting of a static universality class into two dynamical ones 
[92, 17] and, for the case of crosslinking long polymer-chain molecules (vulcanisation), a 
decrease of the width of the critical region with increasing chain length [108]. The latter may 
explain the observation of a crossover behaviour to mean-field properties in certain gelation 
experiments, if measurements were not performed well inside the true critical region. 

As far as we know, no measurements of the critical behaviour of the first normal stress 
coefficient have been reported. The Rouse value k « 0.71 for the viscosity and three- 
dimensional percolation statistics agrees well with the experiments of [79, 107, 80, 83] (only 
[83] was included in Table 3 to demonstrate the broad scatter of the viscosity data). On the 
other hand, it is not compatible with the possibly oversimplifying albeit attractive proposal 
[92, 17] to interpret the wide variation of the viscosity exponent fc as a signature of a splitting 
of the static universality class of gelation into different dynamic ones. Indeed, Rouse and 
Zimm dynamics are considered [109,48] to be at the extreme ends of the strength of the 
hydrodynamic interaction. Since the Zimm model does not even predict a divergence at 
the transition, the actual value of k should then lie below the Rouse value according to that 
proposal. 

Hence, the broad scatter of the experimental data calls for additional relevant interactions 
beyond those accounted for in the Zimm or Rouse model. This may be due to the preaveraging 
approximation. In particular, it throws away hydrodynamic interactions among different 
clusters. But we do not expect this to be the sole relevant simplification of the Zimm model, 
because hnear polymers show a decrease in the viscosity when abandoning the preaveraging 
approximation [110], and effects of preaveraging for branched molecules are even more 

§ Unfortunately, the value of b,, resulting from this scaUng argument in the case of ErdSs-R^nyi random graphs 
was incorrectly ascribed to <i = 6 dimensions in the second last paragraph of [42], leading to the wrong statement 
bn = 1/2 there. 
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pronounced than those for Unear ones [111]. Rather it seems that there are no satisfactory 
explanations without considering excluded-volume interactions. Indeed, simulations [28] of 
the bond-fluctuation model deliver higher values k ~ 1.3 in accordance with the scahng 
relation k = 2i/—f3, which arises from heuristically merging Rouse-type and excluded-volume 
properties, see above. On the other hand, entanglement effects are neglected, too. These 
topological interactions are argued to play a vital role in stress relaxation. However, temporary 
entanglements are expected to play only a minor role [84] for the dynamics close to the 
gelation transition. This is because the time scale of a temporary entanglement is determined 
by the smaller clusters, whereas near-critical dynamics is determined by the largest clusters, 
which contribute the longest time scales. Yet, there remain permanent entanglements due 
to interlocking loops. They are clearly far beyond the scope of the present and many other 
theoretical approaches. 

5. Closing remarks 

The list of shortcomings of the Rouse and the Zimm model for crosslinked monomers is 
long, and it was discussed in Sections 3.3 and 4.3. Yet, one should not underestimate the 
importance of these models for our understanding of gelling liquids. First, the success of 
Theoretical Physics and, in particular. Statistical Physics has always relied on capturing the 
essence of observable phenomena in simple mathematical models. Models that isolate certain 
physical mechanisms and, at the same time, sacrifice many details of the observed reality. It 
is safe to say that, at least for linear polymers, the Rouse and the Zimm model have proven 
to be among this class [47,48]. Second, simple exactly solvable models always represent 
cornerstones against which more elaborate theories, approximation methods and numerical 
simulations can be tested. Moreover, in the absence of an ultimate theoretical picture, the 
predictions of such minimal models also serve as a standard reference for experimental data. 
Indeed, this has been common practice in experimental investigations on gelling liquids over 
the years, see e.g. the review articles [1 12, 44]. All the more it is important to have reliable 
and mathematically firmly based predictions of these model. 
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